Finite temperature phase diagram of ultrathin magnetic films without external fields 
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We analyze the finite temperature phase diagram of ultrathin magnetic films by introducing a 
mean field theory, valid in the low anisotropy regime, i.e., close to de Spin Reorientation Transition. 
The theoretical results are compared with Monte Carlo simulations carried out on a microscopic 
Heisenberg model. Connections between the finite temperature behavior and the ground state 
properties of the system are established. Several properties of the stripes pattern, such as the 
presence of canted states, the stripes width variation phenomenon and the associated magnetization 
profiles are also analyzed. 
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I. INTRODUCTION 

Despite the increasing growth of knowledge about mag- 
netic ordering in ultrathin magnetic films during the 
last decade, both from experimental^— and theoretical 
work a 7 i n ~ — , there are still many open questions, spe- 
cially regarding its finite temperature behavior. One of 
the main obstacles to advance in these studies is the long 
range character of the dipolar interactions, which are fun- 
damental to explain pattern formation in those systems. 
In particular numerical simulations, although have been 
of great aicU^r— , are strongly limited by finite size ef- 
fects. To avoid them, system size must be large enough 
to contain a large number of domains. The main prob- 
lem relies not on the direct influence of dipolar interac- 
tions on the boundary conditions, but on the fact that 
the basic spatial scale for these systems, namely the typ- 
ical domain size, scales exponentially with the exchange 
to dipolar couplings ratio 8 at very low temperatures 2 ^, 
and roughly linear with 6 close to the transition to a dis- 
ordered stat o 2 i 7 i 27 . Typical values of 8 in ultrathin mag- 
netic films, like Fe based films, are around^ 8 ~ 100, thus 
implying the necessity of very large system sizes to acom- 
modate a reasonable number of domains. To perform 
simulations with those sizes represents up to now, even in 
the best case (close to the transition) , a formidable task. 
Therefore, knowledge about how the different thermo- 
dynamical properties scale with 5 would be very helpful 
to estimate whether the numerical results for relatively 
small values of 8 (typically between 3 and 5 up to now) 
can be extrapolated to more realistic values. 

For the analysis of the magnetic properties of ultrathin 
films, the out of plane anisotropy to dipolar coupling 77 
is also important. The system behaviour appears to be 
strongly dependent on experimental features that modify 
it, such as the film thickness and the sample preparation 
conditions. A strong dependence is also observed in nu- 
merical simulations for— 8 = 3. For low r/ values, a SRT 



from a uniformly magnetized planar phase into a perpen- 
dicular striped phase can happen at finite temperature T 
(in the absence of an external field) , in agreement with a 
previous theoretical prediction 2 ^. On the other hand, for 
high values of 77 there is no planar ferromagnetic phase 
and the system undergoes a direct transition from the 
striped state into a disordered one. From these numer- 
ical results a global (77, T) phase diagram was obtained, 
that is in qualitative agreement with a variety of exper- 
imental results 2 ^. However, for such a small value of 8 
certain features can be very different from that expected 
for large values of 8. For instance, at zero temperature 
the striped equilibrium state for 77 above certain critical 
value r\ c (where the SRT occurs) is characterized by a 
stripe width almost independent of 77 for 8 < 5. On the 
contrary, for values of 8 > 5, a strong variation of the 
equilibrium stripe width with 77 emerges wheni^ 77 > i] c . 
Another feature that depends strongly on the interplay 
between exchange and anisotropy is the structure of the 
magnetization pattern close to the SRT. At zero temper- 
ature and close to the SRT, the out of plane component 
of the magnetization presents an almost sinusoidal shape 
with a large in-plane component, displaying a canted 
structure^. For small values of 8 (5 < 5) such struc- 
ture remains for a relatively large interval of values of 77 
above the SRT and changes abruptly to a completely per- 
pendicular striped state with sharp domain walls (Ising 
like state). Consistently, numerical evidences of a canted 
structure with a sinusoidally shaped magnetization pro- 
file close to the SRT at finite temperature has been re- 
cently reported 2 ^ for 8 — 4.5. However, as 8 increases the 
range of anisotropy values at which such canted state is 
present shrinks at zero temperature^, becoming almost 
negligible for realistic values of 8. Hence, it is not clear 
whether it is expected to be relevant at finite temperature 
or not. 

In this work we analyze the finite temperature phase 
diagram in the low anisotropy region (close to the SRT) 
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and several related properties using a coarse-grained 
based mean field model for ultrathin magnetic films and 
Monte Carlo simulations on a microscopic model. The 
main objective of the paper is to discuss which of the ob- 
served features of the phase diagram for low values of S 
are expected to reflect the large S behavior. Several prop- 
erties stripes patterns are also analyzed. The plan of the 
paper is as follows: in Section [II] we introduce the coarse 
grained model and calculate the associated mean field 
phase diagram. In Section IIIII we present Monte Carlo 
simulations results for a Heisenberg model and compare 
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them with the previous ones. In Section IIVI we discuss 
our results. 



II. THE MEAN FIELD MODEL 

We consider a general phenomenological Landau 
Ginzburg free energy for a two dimension ultrathin mag- 
netic film of the form, 



F[M] 



X - J d 2 x{(VM(x)) 2 +r M 2 (x) + ^M 4 (x)} 
-|jd 2 xM 2 2 (x), 
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d 2 x / d 2 x' 



M(x).M(x') - 3(n.M(x))(n.M(x')) 



(1) 
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where M = (M x , M y , M z ) is the coarse grained magneti- 
zation. 6 is the exchange to dipolar couplings ratio, rj is 
the anisotropy to dipolar coupling ratio and n is a unit 
vector pointing in the x — x' direction. A cutoff at some 
microscopic scale A is implied in the second integral. We 
will assume A = 1. The temperature dependency comes 
through ro = ro(T). 

In order to minimize Eq.([T]), we propose a variational 
stripe like solution, i.e., a modulated solution along the 
y direction, where only Bloch walls between domains are 
allowed 2 ^, namely M(x) = M(x) and M x (x) = 0. We 
also assume that modulus of the magnetization is uni- 
form, i.e. 



M 2 (x) + M 2 {x) = M 2 Vx. 

This approximation is expected to breakdown for large 
enough values of rj, where the statistical weight of spin 
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configurations with large in plane components tends to 
zero, but non uniform out of plane configurations are still 
expected to minimize the free energy^. In fact, in the 
77 — > 00 limit the whole effective free energy ([1]) cease to 
be valid, being replaced by a functional of a scalar order 
parameter (local out of plane magnetization), without 
the anisotropy terrr*i£. 

Under the present assumptions, the following form of 
the dipolar term can be assumed 2 ^ 



1 1 



dx 



,M z {x)M z (x>) 



where we have neglected the self energy term arising from 
the dipolar energy, since it just implies a constant shift 
in the anisotropy coefficient rj. 

Then, the variational free energy per unit area reduces 
to 



,r,xi 1 f , \fdM v \ 2 fdM z \ 2 ,,, u .1 I f 1 f , .Mjx)MJx') k 

m = uj*\{i*r) + UrJ + roM + 2 \ + iLj dx J dx l - 4 l*. 1 dx (2) 
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wher&ili k = rj — a with a = 3.485 . . . We can write where, 
M z (x) = M<p(x) where \(j){x)\ < 1. Then 



f=l(r (T)+2e/S) M 2 + ^M* (3) 
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dx l dx 



(x — x') 2 



dx (j) 2 (x), 



(4) 



i.e. e[(/)(x)] is the microscopic energy per spin of the 
associated microscopic model (Heisenberg model with 
out of plane anisotropy, exchange and dipolar interac- 
tions), for a spin density profile (S x (x), S v (x), S z (x)) = 
(0, yl — 4>{x), Its minimal energy configuration 

as a function of the microscopic parameters (5, 77) can be 
described by variational expressions characterized by dif- 
ferent sets of variational parameters Pi,P2, ■ ■ ■ that will be 
described later. Minimization of the free energy Eq.Q 
leads to: 



df 

= M(r +2e( Pl , P2 ,...)/S + uM 2 ) 



df_ 

dpi 

We also have; 



M 2 de 
8 dpi 



d 2 f 
dM 2 



= r + 2e/S + 3uM 2 



= (5) 
(G) 

(7) 



and 



M 2 = -{r + 2e/8)/u, 



de 
d Pj 



0. 



(11) 
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Hence, from Eqs.© and pT|) it follows that, 



d 2 f 
dM 2 



2uM z 



and from Eqs.© and (fT2|) . 

d 2 f 



dMdpi 



0. 



and 



d 2 f 
dMdpi 
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5 dpi 

d 2 e 



-SHI— = I ]\r 2 _ 

dp l dp ] 8 dpidpj 



(8) 
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M = is always a solution of the extremal equations 
(H])-© and the corresponding free energy / = is inde- 
pendent of the parameter values of (pi,p2, ■ ■ .). Hence, 
all the second derivatives are zero, except 



d 2 f 



dM 2 



= r 



2e/S, 



(10) 



M=0 



which controls the stability of the M = solution. 

An ordered solution (local minimum of /) with M =/= 
exists whenever ro + 2e/8 < 0. From Eqs.([5])-([6]) we have 
that: 



Then, the Hessian matrix of / has a positive eigenvalue 
2u M 2 and a diagonal block that equals M 2 H e /8, where 
H e is the Hessian matrix of e. Therefore, a local min- 
imum of / has to be a local minimum of e. The free 
energy of an ordered phase is: 



f = -—(r a (T) + 2e mm /8) 2 . 

4:11 



(13) 



The extremal properties of e are well know n 16 ' 27 . For 
low values of the anisotropy rj the minimum of e cor- 
responds to a planar ferromagnetic (PF) configuration 
4>(x) = 0. Above certain critical valu o 16 ' 27 rj c (8) = 
a + 7r 2 /3 — 7r 2 /2<5, e is minimized by a striped profile with 
periodicity 2h (h is the stripe width). Close to r\ c the do- 
main structure corresponds to a canted sinusoidal wall 
profile (SWP), where \<\>{x)\ = cos (9 is constant inside 
the striped domains (0 < < ir/2 is the canting angle) 
and presents a sinusoidal structured wall of width^l w. 
The energy of the SWP is given by 
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e SWP ( S , k, A) = ^ (l - VT^ 2 ) + s 2 (V /3 - V) (1 - A/2) - ^ G(A)) , 



(14) 



where s = cos 6*, k = n/h, a A = w/h, ands 7 - 
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G ( A ) = "a Z2 —R 

7r z * — ' m ( 1 — r, 



m— 1,3,... 



7rmA 

{l-m 2 A 2 ) 2 C ° S " \~ 



(15) 
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FIG. 1: Mean field phase diagram for S — 6 and cl/Tf = 1. 
Solid lines correspond to second order phase transitions. The 
dashed line marks the crossover between canted stripes and 
perpendicular stripes configurations; it is estimated arbitrar- 
ily as the region above which the maximum in-plane compo- 
nent of the magnetization is less than 5% of the saturation 
magnetization M (s < 0.05) 



Different approximations of high accuracy for G(A) are 
available^^, so the values of (s, k, A) that minimize 
Eq.(|14p can be found numerically for arbitrary values of 
(5, rj) . As the value of r) is raised above r\ c the canting 
angle decreases fast from 9 = ir/2 at rj = r\ c to 9 m and 
the striped pattern that minimizes e changes to a hyper- 
bolic wall magnetization profile (HPW) whose energy is 
given byi£ 

e HPW (k,A)= 1 (l-A/2) + — —-— In^J , (16) 

where 7 = A - 77, with A = 4.5327.... Ed.([T6]) can be 
easily minimized^ 

If r < and 77 < r\ c (e min = 0), the global minimum 
of / corresponds to M 2 = — ro/u and s = 0, that is, 
to a planar ferromagnetic (PF) state with free energy 
/ = — Tq/4u. When ro = and rj < ?/ c the system 
undergoes a second order phase transition between the 
paramagnetic state and the PF one, independently of 77. 
We will assume hereafter that ro = a(T — Tf) , where Tf 
is the paramagnetic to PF transition temperature. 

When T < Tp and 77 > ?7 C the SWP configuration with 
free energy given by Eq. (|13|) is the stable solution for 



values of rj close to ?7 C . Since the striped order emerges 
continuously, the SRT at 77 = r/ c , according to the present 
approximation, is a second order one. As 77 is further in- 
creased there is an energy crossing at certain value of 
77 and the stable configuration changes into an HWP. 
Hence, for rj c < < rj* the stable configuration is a 
canted striped one, while for 77 > 77* the stripes are fully 
saturated in the out of plane direction inside the domains 
(we will call this state an "Ising striped configuration"). 

If T > Tp (ro > 0) and rj > rj c the global minimum cor- 
responds to the modulated phase when ro + 2e TO j n /S < 0, 
i.e. T > Tp — 2e m i n /aS. Therefore, there is a transition 
line at T c (rj) = Tf — 2e m in{rj) / a8 . The order parameter 
changes continuously at T c (M 2 = — (ro + 2e TO , n ) = 0), 
but s changes discontinuously. In Fig[T]we illustrate the 
typical topology of the phase diagram for the particular 
case 8 = 6. All the solid lines in Q] correspond to sec- 
ond order phase transitions. We also show the crossover 
line between the region where the magnetization profile 
shows a significative canting angle (canted stripes) and 
the region where the local magnetization is almost per- 
pendicular to the plane (Perpendicular stripes). 

Although the paramagnetic solution is not a global 
minimum of / when rj > rj c and T > Tp, it could be 
still a local minimum provided that r<j + G < for some 
values of s and ko ■ From Eq. (fT0"|) , such condition ensures 
the local stability against variations of M. However, since 
all the rest of the second derivatives cancel, the complete 
stability of the paramagnetic solution is beyond the linear 
analysis. We verified numerically that indeed the param- 
agnetic phase remains locally stable (metastable) below 
T = Tp at fixed rj > rj c . Such metastability is a result of 
the high degeneracy of the paramagnetic solution under 
the present approximation, so it appears to be a spurious 
result. However, it can be indicative of a change in the 
order of the transition if the approximation is improved. 
Indeed, there are several evidences towards the first order 
nature of the stripes-disordered phase transitio n 11 ! 19 ' 23 . 



III. MONTE CARLO SIMULATIONS 

In order to compare the mean field results with the 
behavior of a specific microscopic model, we performed 
Monte Carlo (MC) simulations using a Heisenberg model 
with exchange and dipolar interactions, as well as uni- 
axial out of plane anisotropy. The model, which de- 
scribes an ultrathin magnetic film (see Refll6l and refer- 
ences therein) can be characterized by the dimensionless 
Hamiltonian: 



Sj, • Sq 



- 3 



[Si ■ fij) {Sj ■ fij) 



(17) 



<»,i> 



where the exchange and anisotropy constants are normal- 
ized relative to the dipolar coupling constant, < i,j > 
stands for a sum over nearest neighbors pairs of sites in a 
square lattice with N = L x x L y sites (the lattice param- 
eter is taken equal to one), stands for a sum over all 



distinct pairs and 



is the distance between 



spins i and j . Each spin is defined by a unit vector with 
components S X ,S V ,S Z . All the simulations were done 
using the Metropolis algorithm, and periodic boundary 
conditions were imposed on the lattice by means of the 
Ewald sums technique. We focus our simulations on the 
case 6 = 6, where the system presents a canted equilib- 
rium state at zero temperature for a wide range of the 
anisotropy values^. 

The phase diagram was obtained by measuring the out- 
of-plane magnetization; 



M, 



(18) 



(the in plane components are defined in a similar way) 
the in-plane magnetization; 



and an orientational order parameter—; 

n h - n v 



Ohv 



(19) 



(20) 



where (• • • ) stands for a thermal average, (n„) is the 
number of horizontal (vertical) pairs of nearest neighbor 
spins with antialigned perpendicular component, i.e., 



= ^E {l-sig[S z (r x ,r y ),S z (r x + l,r y )}} (21) 



and a similar definition for n v , where sig(x, y) is the sign 
of the product of x and y. To obtain the stripe width of 
the modulated states we considered the structure factor 
\S{k)\ 2 , where 



S(k) 



1 



-ik.r 



(22) 



The stripe width was calculated using the expression 



ir/k r , 



(23) 



where k max is the modulus of the wave vector that max- 
imize |S(fc)| 2 . 

In order to find the equilibrium phase diagram, we car- 
ried out the simulations with two protocols for the inde- 
pendent parameter (temperature, anisotropy or external 
field). In the first protocol, we varied the independent 
parameter linearly with the simulation time, increasing 
or decreasing it at a given rate r, keeping the rest of the 
parameters fixed. For instance, if we choose the temper- 
ature, then T(t) = T(0) + rt where t is the simulation 




FIG. 2: (Color online) In plane magnetization and orienta- 
tional order parameter as a function of anisotropy for 5 = 6. 
T = 0.1 (upper panel) and T = 1.0 (lower panel). 



time measured in units of Monte Carlo Steps (MCS). 
Each MCS corresponds to N single spin updates of the 
Metropolis algorithm. The initial spin configuration at 
T(0) was previously obtained by performing t e MCS to 
equilibrate. The order parameters were calculated along 
the simulation and averaged over many realizations to 
improve statistics. We called this protocol "linear varia- 
tion of parameters" (LVP). The second one was a ladder 
protocol. For instance, in the case of T being the inde- 
pendent parameter, the system is initialized at the para- 
magnetic state at a high temperature, and then tempera- 
ture is reduced at discrete steps. The initial configuration 
for each temperature is the last one of the previous step. 
At each step we discarded the first t e MCS in order to 
equilibrate, then we calculated the averages over the next 
t m MCS. 

First, we calculated M\\ and Ohv as a function of the 
anisotropy for two fixed temperatures. We applied the 
LVP protocol to increase r\ from a small value r\ < ?/ c , 
starting from an equilibrated in-plane ferromagnetic con- 
figuration. In these simulations we used L x = L y = L, 
with L = 80 and 120. The parameter variation rate r 
was ranged from r = 10~ 5 to 10~ 7 , depending on the 
temperature and the system size. The typical behavior 
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FIG. 3: (Color online) Orientational order parameter (red cir- 
cles) and in-plane magnetization (black squares) as a function 
of temperature for 5 = 6 and r\ = 8.7. 



of the order parameters at low temperatures (T = 0.1 
and T = 1.0) is illustrated in FigH Three different 
behaviors can be identified: the planar ferromagnetic 
state at low anisotropics, characterized by M\\ ^ and 
Oh v = 0; a canted striped state at intermediate values of 
77, with both M|| ^ and Ohv 7^ 0; and a perpendicular 
striped state at large enough values of 77, characterized 
by M11 = and Ohv 7^ 0. At each temperature the tran- 
sition points are identified with the values at which the 
corresponding order parameter becomes zero. 

At large anisotropy values (77 > 8.4). the in-plane fer- 
romagnetic phase is absent for any value of the temper- 
ature. The system undergoes a direct transition from an 
almost perpendicular striped state into the paramagnetic 
state as the temperature is increased. This can be seen in 
FigJ21 where the order parameters were computed using 
a ladder protocol cooling from T = 4.0, with L = 120, 
t e = 10 5 , and t m = 10 5 . 

By means of the methods described above, we obtained 
the phase diagram shown in FigJU Although, the global 
topology of this diagram is in agreement with that ob- 
tained by the mean field theory (see Figfl]), some notice- 
able differences exist, which appear to be an artifact of 
the mean field approach. Some of them will be discussed 
111 section HVl The presence of a canted striped region is 
in agreement with previous MC calculations carried out 
by Whitehead et al. in Ref. [24| for 6 = 4.5 and with the 
zero temperature behavior of the models. In particular, 
our results show that the canted region is larger than the 
corresponding to 6 = 4.5. 

We explored the canted region looking at the variation 
of the stripes width h. This is a difficult task because 
the stripe width variation is mediated by the formation 
of topological defects (usually stripe dislocations) which 
need long simulation times to nucleate and move. An 
acceleration of this process was obtained when we added 
an in-plane external field term to Eq. (|17[) of the form 
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FIG. 4: Monte Carlo phase diagram for 5 = 6. 



— £$^(<Sf) an d applied a LPV protocol with £ as the 
free parameter. The value of £ vary from £(0) = H x 
where H x is an external magnetic field strong enough 
to saturate the magnetization in the x direction, to zero 
(zero-field condition). After several tests, we found that 
H x = 0.5 was optimal for all the regions of the phase 
diagram studied in this work. Then, tu extra MCS were 
performed before calculating h. The stripe width shown 
as a function of temperature for 77 = 7.5 in FigJ5]is the 
result of an average performed over several realizations 
of the LPV protocol and the error bars correspond to the 
dispersion of h. 

The stripe width at the SRT is h = 10 and remains 
constant down to T = 0.75, where it displays a steep 
increase. The widest stripes are observed at low temper- 
atures reaching a maximum value of h ~ 14, close to the 
zero temperature value, h = 17, calculated previously^. 
The h values obtained for both system sizes are almost 
undistinguishable, showing that finite size effects are neg- 
ligible. 

We next analyzed the stripe width variation with the 
anisotropy, which is closely related to the variation with 
the film thickness d. Indeed, previous numerical simula- 
tions suggest that the effective out of plane anisotropy 
varies inversely with the film thichness^ 77 ~ 1/d. A 
usual experimental technique to analyze the effects of the 
film thickness on the magnetic patterns is to take images 
on wedge-like ultrathin films, where the film width is a 
function one of the axis (eg. RefsilH). We modeled 
these systems assuming the relation 77 = l/d(x), where 
d(x) is the local film width that depends on the x posi- 
tion. The simulations were performed over rectangular 
lattices of size L x = 360, L y = 180 using the external 
field protocol used in the study of the stripe width varia- 
tion and periodic boundary conditions in the y direction. 
The results are shown in Fig HI] in two columns, each col- 
umn being related to different d(x) assumptions. The 
left one corresponds to a continuous linear variation of 
d(x) = a + bx, where a and b were chosen such that 77 
varies from 77 = 6 at x = 1 to 77 = 8 at x = 300. The 
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FIG. 5: (Color online) Stripe width as a function of temper- 
ature for 5 = 6 and r\ = 7.5. The simulation parameters were 
L = 80, r = -10" 6 , i e = 10 5 and t H = 10 5 (circle) and 
L = 120, r- = 10 -7 , t e = 10 6 and t H = 10 6 (box). The er- 
ror bars are the standard deviation taken into account many 
realizations, some of them are omitted for clarity. 



right column corresponds to a ladder structure, where 
d(x) varies at equal-spaced steps corresponding to the 
values rj = 8, 7,6 from left to right. The figure presents 
typical snapshots of equilibrated magnetic patterns at 
fixed temperature. Each pixel represents a spin compo- 
nent in gray scale, ranging from white when the value is 1 
to black when it is — 1. From the S z component behavior 
turns out that the stripes width decreases as x increases 
until the SRT. Once the SRT is reached, the spins are 
ferromagnetically ordered in the same direction of the 
in-plane component of the magnetization in the walls. 
Moreover, the S x components in the walls are along the 
stripes direction, showing the they are Bloch's walls as 
expected^. The stripe width reduction occurs by the 
insertion of new stripes from the low anisotropy (higher 
thickness) region, in agreement with experimental results 
on Fe on Cui and Fe/Ni on Cu3 films. These results give 
further support to the assumption r) ~ 1/d and suggest 
that the observed stripe width variation with the film 
thickness is due to the induced anisotropy gradient. In 
fact, the r) ~ 1/d dependency is probably related to the 
contribution to the effective anisotropy coming from the 
short range part of the dipolar energy, which can be as- 
sumed proportional to the film thickness^ (at least in 
the ultrathin limit). The stripe width variation with rj 
in the wedge like film of Fig. [5] is shown in Fig. [7] We 
also performed a series of simulations on lattices with 
L x = L y and uniform anisotropy, for different values of 
rj. The equilibrium average stripe width agreed with that 
observed in the wedges. 

We verified that the increase in the stripe width as 
r] increases follows a series of steps in a similar way to 
that observed at zero temperature^. In other words, the 
canted region in the phase diagram of Figgis composed 
by a series of transition lines (not shown for clarity) that 
follow a similar direction as the SRT line and converge to 





FIG. 6: Snapshots of the spins components for two wedges in 
equilibrium for 8 = 6 at T = 1. The film width functional 
dependence varies as a linear function of x (left) or as equal- 
spaced steps of values rj — 8,7,6 (right), as schematized at 
the bottom of the columns. The simulation parameters are 
r = 10" 5 and t e = I0 5 . 



the zero temperature transition points between different 
stripe width ground states 1 ^. Every time the anisotropy 
crosses one of such lines the stripe width increases by 
one unit (dynamically mediated by defects). In this way, 
the stripe width variation with temperature (horizontally 
crossing of those lines in Fig|2} is related to the ground 
state structure of the system. In this context, the absence 
of stripe width variation with temperature for rj > 8.5 is 
related to the fact that for 8 = 6 the ground state stripe 
width has already saturated 1 ^ 

Finally, we analyzed the magnetization profile varia- 
tion of the stripe pattern as a function of anisotropy, 
which can be seen in Fig GO These results were obtained 
by averaging over 25 adjacent profiles S z (x,yo) in a sys- 
tem with uniform anisotropy in equilibrium at T = 0.5. 
Topological defects like dislocations were avoided in the 
calculation. We thermalized the system using the ex- 
ternal field protocol with parameters L x = L y = 144, 
r = 1(T 5 , t e = 10 5 and t m = 10 4 . At 77 = 6.45 the 
system is in the canted state, the walls are wide with 
sinusoidal like shape. The perpendicular components of 
the spins at the center of the stripes is lower than 1, 
meaning that all the spins have an in-plane component 
aligned with the stripes, thus contributing to the in-plane 
magnetization as pointed out in Figj2j As the anisotropy 
increases, the spins at the center of the stripes become 
perpendicular to the plane and the walls narrow but are 
still extended. The same behavior of the magnetization 
profile is observed at zero temperature 1 ^. Finally, at high 
enough anisotropy values (77 > 8.0), the stripes widen and 
the wall widths w become close to one. 
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FIG. 7: Stripe width as a function of 77 corresponding to the 
wedge like film of the left panel of Fig(6] (5 = 6). The dashed 
line correspond to a parabolic fiting. 




FIG. 8: (Color online) Mean stripe magnetization profiles in 
the perpendicular direction for S — 6 at T — 0.5. Each curve 
is an average of consecutive profiles along a stripe and the 
error bars correspond to the standard deviation. 



It is worth to note that the fluctuations of the spin 
directions at low and intermediate anisotropy values are 
stronger within the walls, as can be observed from the 
error bars. This suggests that these spins are less re- 
stricted to move and therefore facilitate defects mobility. 
This explains the higher efficiency of the previously used 
in-plane field protocol to obtain thermal equilibration, 
through the interaction between the external field and 
the large in plane components inside the walls. 

Interestingly, the change in the magnetization stripe 
profile as the anisotropy increases closely resembles 
that observed as the temperature decreases, both 
experimentally 7 - and in mean field theorie a 7 ' 17 . 



IV. DISCUSSION 

Our mean field results suggest that the global topol- 
ogy of the phase diagram observed numerically for low 



values of 5 (both from previou o 23 ' 24 and from the present 
simulations) is robust, at least under the validity condi- 
tions of the present approximation, namely, in the low 
anisotropy region close to the SRT. For large enough val- 
ues of the anisotropy the approximation breaks down, as 
evidenced by the unphysical monotonous increase of the 
transition temperature between the stripes and param- 
agnetic phases in the large 77 region. This breakdown is 
on the basis of the present MF approximation, namely, 
in the effective free energy Eq.([T|). Such free energy 
can be obtained variationally from a partition function 
Z = Tre~^ H [ M l, where the coarse-grained Hamiltonian 
if[M] has the same structure of 3 - 1 - Eq.lJTJ. The effective 
free energy -F[M] is then the order zero term in an expan- 
sion of i?[M] around its minimum when fluctuations are 
neglected. If the 77 — ¥ 00 limit is taken a priori of such 
expansion, all the configurations with non zero in plane 
magnetization components get zero statistical weight and 
a different Landau Ginzburg free energy (which depends 
only on the scalar field M z ) is obtained^ 7 -. Therefore, 
the correct stripe-paramagnet critical temperature must 
converge to the (77-independent) value predicted by the 
last free energy when r\ 3> 1. In other words, even within 
the mean field theory the correct behavior cannot be ob- 
tained as the 77 — > 00 limit of the present approach. 

While the previous difference (vertical line vs. finite 
slope) between the MC and the MF phase diagrams is 
particular of the present approach, some others appear to 
be associated to general features of the mean field theory. 
For instance, the transition line between the planar fer- 
romagnet and the canted-stripe phases computed within 
the MF approximation is horizontal, while it shows a fi- 
nite slope when extracted from MC simulations. This 
seems to be a direct effect of neglecting thermal fluctua- 
tions, since theoretical works show that those fluctuations 
renormalize the dipolar and anisotropy coupling param- 
eters in such a way that the anisotropy K(T) diminishes 
faster than the dipolar coupling constan t 28 ' 32 g(T) (in 
our notation, 77 = K/g). Those works predict a linear 
dependence of the reorientation transition temperature 
with anisotropy with positive slope, which is roughly in 
agreement with the transition lines obtained from MC 
simulations. Finally, the transition line between the pla- 
nar ferromagnet and the paramagnetic phases is a verti- 
cal straight line in MF diagram while it shows some slope 
in the MC diagram. This is because of the simplifying 
assumption that the (coarse grained) phenomenological 
transition temperature Tp is independent of 77. While the 
tendency of the transition line in the MC diagram sug- 
gests that this assumption may appropriately describe 
the very low r\ limit, it clearly fails for large enough val- 
ues of 77. An increase in the perpendicular anisotropy 
should destabilize the planar ferromagnetic phase, thus 
decreasing the critical temperature. 

One fact that emerges, both from our mean field and 
Monte Carlo results, is the strong influence of the ground 
state properties on the finite temperature behavior close 
to the SRT. One example is the presence of canted states 
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close to the SRT line. Comparing with previous MC 
results for— 8 = 4.5, the phase diagram canted region 
becomes wider for 5 = 6, consistently with the zero tem- 
perature phase diagram^. However, the range of values 
of r\ where the ground state canted angle is different from 
zero becomes extremely narrow as 5 is further increased. 
Hence, our mean field results suggest that those states 
would be present at finite temperature only very close 
to the SRT line for any realistic value of 5. Another 
example is the stripe width variation and the magnetiza- 
tion stripes profile change with 77, that closely follow the 
zero temperature behavior—. Moreover, the qualitative 
agreement between our MC simulations on wedges and 
experimental results on Fe/Ni ultrathin films^ supports 
inverse relationship between out of plane anisotropy and 
film thickness rj ~ 1/d. 

The correlation between the stripe width variation 
with the temperature and with the anisotropy observed 



close to the SRT in the present simulations is another 
interesting fact. As previously pointed oufei, varying the 
film thickness (always in the ultrathin limit) produces 
a similar effect as changing the temperature, thus lead- 
ing to an "inverse effective temperature" interpretation 
of the thickness^. Considering the relation 77 ~ 1/d, this 
appears to be consistent with the similarity observed be- 
tween the change in the magnetization profile when 77 is 
varied and that observed in Fe films when the temper- 
ature is varied^. Such set of similarities suggest that a 
deeper analysis about the interplay between temperature 
and anisotropy could shed an additional light about the 
origin of the strong stripe width variation with tempera- 
ture observed in ultrathin magnetic films. 
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